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I give a review and progress report on studies of lattice QED. I emphasize analytical results and methods that 
are applied in data analysis. Also, I derive some bounds for the critical exponents and establish their connection 
with scaling violations. Triviality, as realized in (f> 4 theory, is ruled out on theoretical grounds. I show that the 
present data, if analyzed correctly, all lead to the same conclusions. They are compatible with power law scaling 
with nongaussian exponents. 



1. Introduction 

The most important theoretical issue related to 
QED is the question of its existence and the na- 
ture of the continuum limit. Traditional wisdom 
is that non-asymptotically theories do not have 
an interacting continuum limit in four dimension. 
This belief has been based on perturbation theory 
and, until recently, hasn't been challenged outside 
of the pcrturbative context. The study of lattice 
QED is an atcmpt to clarify this issue Q || || Q . 

The central problem in the study of the con- 
tinuum limit of a theory is the verification of hy- 
perscaling ||. The starting point is the free en- 
ergy from which all connected correlation func- 



tions can be generated: F s -, 



= i A 



*F(h/t*). 



Having the dimension of inverse volume, hyper- 
scaling implies that F s i ng ~ £ . The renormal- 
ized coupling is defined as 
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where x ( ™° = d 3 M/dh 3 is the zero- momentum 
projection of the connected four-point function. 
If the action is gaussian, Wick's theorem applies 
and the iV-point functions factorizc. Thus, for 
a gaussian theory g^ vanishes and is non-zero 
otherwise. Using the hyperscaling hypothesis, 
eq.(l.l) can be converted into 



(2A-7— dv)/v 



(1.2) 



Being dimensionless, should be independent 
of £ if £ is the only scale. Thus, the validity 



of hyperscaling requires that the exponent must 
vanish. So, hyperscaling implies the relation, 
2A — 7 — dv = 0, between the critical indices. 
It is known in general that the following inequal- 
ity holds 2 A < 7 + dv. This means that the 
exponent in the expression for gn is always non- 
positive, so that violations of hyperscaling im- 
ply that the resulting theory is non-interacting. 
Above four dimensions, the exponents are gaus- 
sian (7 = 1, A = 3/2, v — 1/2). In this case, it is 
easy to verify the above inequality - it amounts to 
d > 4. In four dimensions most field theoretical 
models critical exponents have the mean-field val- 
ues, but with logarithmic corections that drive ga 
to zero. From here, the importance of the knowl- 
edge of critical exponents becomes clear. Estab- 
lishing their non-gaussian values would be a sig- 
nificant step towads ruling out the triviality of 
the theory. 

The recent discovery that non-compact QED 
undergoes a second order phase transition (|] 
from a massless to a massive phase with spon- 
taneously broken chiral symmetry reopened the 
old question of the existence of quantum electro- 
dynamics. In what follows I will survey what has 
been done. 

When studying chiral symmetry breaking in 
gauge theories one is faced with the following 
problem. Because of the singular nature of the 
chiral condensate, the chiral limit is not directly 
accessible in numerical simulations and work at 
finite bare mass is required. In this way the pre- 
cise position of the critical coupling is difficult to 
determine. 
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At the critical point eq.(1.4) results in a simple 
equation 
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(1.5) 



The ratio is independent of the symmetry break- 
ing field! The knowledge of the mass ratio can 
be used to determine both the critical point and 
exponent S B. 



Figure 1. Mass ratio versus bare mass for N = 2 
theory on a 16 4 lattice. Solid lines are the fits 
obtained from the EOS. 



A usual way to procede is to commit oneself to 
some scenario and extrapolate to the chiral limit 
looking for the best fit etc.. The index S deter- 
mines the approach to the chiral limit. It is de- 
fined at the critical point by m ~< ijjip > s , (t = 
0). Here, I will discuss an alternative to the tradi- 
tional approach to finding the critical exponents 
in a given theory and show how it is applied to 
QED where the traditional approach led to am- 
biguous and controversial results in the past. 

Consider the effect of chiral symmetry break- 
ing on the meson spectrum from a physical point 
of view. In the symmetric phase, chiral symme- 
try requires the degeneracy between the chiral 
partners. Therefore, in the chiral limit the ra- 
tio R — Ml I Ml — 1. As the bare mass increases, 
the ratio decreases. In the broken phase, however, 
the ratio vanishes in the chiral limit because the 
pion is a Goldstone boson. Behavior of the mass 
ratio is displayed in Fig.l. The data are taken 
from ||. 

The connection between the mass ratio and the 
order parameter is obtained from the equation of 
state (EOS). 



h a = M a M 5 - 1 f( 
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2. Two and Four Flavor Theory 

I first discuss the four flavor data. Simulations 
in this case were done on smaller lattices (12 4 by 
DESY § and 10 4 by Illinois group @). For 
N = 2 theory simulations were done on a bigger 
lattice (16 4 ) with better statistics by the Illinois 
group ||. The quality of the four flavor data 
is much more modest and are contaminated by 
the finite size effects. This being the case, the 
conclusions that have been drawn in the past are 
to be taken with caution. 



Figure 2. Mass ratio data versus m taken from 
the DESY data. The simbols are explained in the 
text. 



The plot of mass ratios versus the bare mass 
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is shown in Fig. 2. The data are taken from 
the DESY group 0. Critical parameters from 
these data are extracted simply. They are p c = 
0.205, S — 2.3. The order parameter data, are 
shown in Fig. 3. The data on this plot are from 
the Illinois group [^0[ . They are consistent with 
those of the DESY group, m = 0.03 was the low- 
est mass that was insensitive to the finiteness of 
the lattice Thus, unlike in the analysis of 

the DESY group, m — 0.01 and 0.02 are not in- 
cluded here. Both, f3 c and 6 obtained this way 
agree with the values extracted from the mass ra- 
tios |l(| . It should be noted that the critical cou- 
pling proposed by the DESY group @, /3 C = 0.185 
(dashed line), lies deeply in the strong coupling 
region. Both mass ratio and Fig. 3 show this very 
clearly. No extrapolations or additional assump- 
tions were made so far. Note that if j3 c were in- 
deed 0.185, as proposed by the DESY group, then 
the corresponding isotherm would have to bend 
upwards. Clearly, such a curvature would corre- 
spond to < tpip > that is too small for a given 
value of m, an effect typical for symmetry break- 
ing in a small volume. In the thermodynamic 
limit, however, this is not possible on theoretical 
grounds as I will discuss later. 



predictions for the mass ratio data and the fit 
is given by the solid lines in Fig. 2. The fit is 
compelling and we conclude that consistency is 
satisfactory. 



Figure 4. EOS fit to power-law scaling for N=4 
data from Illinois gruop. 



Now, I review briefly the atempts of the DESY 
group to fit the same data to the logarithmically 
improved mean field theory. In ref. [Q a logarith- 
mically improved 0(2) sigma model is used for 
fitting purposes and it fits the data very well. The 
EOS reads, 



Figure 3. Order parameter data for the four flavor 
theory. 



The EOS (Fig. 4) results in a straight line fit to 
f(x) and exponent (3 = 0.764. This EOS gives 



ln p 



In cr 



(2.1) 



where a = (iptp), t = T]0(1 - f3/f3 c ) and 0- 1 = 
o + 6\{1 — P/P c ). Choosing specific values for 
the five parameters p c , p, t\, 9 c and 9\, a very 
good fit to the data is found, Fig. 5. The result- 
ing chiral transition occurs at p c — 0.186(1) with 
mean field critical indices built in. Since the fit 
involves five parameters, its significance is cer- 
tainly debatable. Let us subject eq.(2.1) to the 
same test that the power-law scaling fit has just 
passed. In particular, from Sec. 8 of ref.[Q we 
read off a formula for R — M%/M% calculated in 
the logarithmically-improved 0(2) sigma model. 
This is done in Fig. 6. Obviously, the fit fails qual- 
itatively. 
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Figure 5. EOS to MF+log scaling, eq.(2.1). 
DESY data. 



Figure 6. Mass ratio data versus the prediction 
from the MF+log EOS. Desy data 



What is wrong with the MF+log scenario? It 
appears to fit the order parameter very well, but 
it fails with the masses. The answer is simple. 
Logarithmic violations of scaling as displayed in 
eq.(2.1) can never appear in the case of chiral 
phase transition. At the critical point, eq.(2.1) 
predicts 

< 4>ii) > 3 



log(l/ < tpip >) 



(2.2) 



Because of the scaling violations, the RHS van- 
ishes faster then a pure power. So the "effec- 
tive" S is bigger then its mean-field value. In the 
Nambu-Jona-Lasinio model for example critical 
EOS gives, on the other hand, 



m 



< t/jip > 3 log(l/ < tpt/; >) (2.3) 



Unlike for magnets, the log's appear in the nu- 
merator - the RHS vanishes slower then the pure 
power and the "effective" S is smaller then the 
(pure) mean-field value. As shown in p2[ , this 
difference in the position of the logarithm i.e. the 
sign of the scaling violations is generic for the two 
models. It is independent of the approximation 
and follows from the basic differences in physics of 
the two systems. In |l^] it was shown that expo- 
nent 8 satisfies two different bounds, depending 
on whehcr Goldstone bosons are elementary or 
composit. For magnets 8 > 3, while for fcrmions 



5 < 3. These bounds are respected by the log- 
arithmic violation of scaling in four dimensions 
p2|| . The origin of these bounds on 8 lies in the 
fact that for magnets weak coupling phase is bro- 
ken wheras for fermions the situation is reversed. 
Thus, the failure of the fits of ref.0 is under- 
standable and expected since it is borrowed from 
the magnetic context and applied to the fermions. 
The source of errors in the fits of ref.|Q] is missi- 
dentifiacion of the physics and the lack of control 
of finite size effects. Because of the small (12 4 ) 
lattice, (3 — 0.185 that was mistakenly identified 
with the critical coupling by the DESY group 
was deep at the strong coupling phase. Such a 
misidentification coresponds to a smaller size of 
the strong coupling phase and the larger value of 
8. 

Concerning the other results reported by the 
DESY group in [Q, several comments are in or- 
der. Some preliminary attempts to extract renor- 
malized charge and fermion mass have been made 
in ref.Q. I should mention some obvious techni- 
cal problems that concern the state of affairs re- 
lated to this attempt. Like order parameter stud- 
ies, simulations of the renormalized charge and 
fermion mass were done on a 12 4 lattice. The pro- 
cedure that was employed assumed the massless- 
ness of the photon. The renormalized charge was 
then extracted as the residue of the photon prop- 
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agator evaluated at zero momentum. However, 
due to the finite volume constraint, the propa- 
gator was evaluated at finite momentum and the 
result was extrapolated to k = 0. This is a myste- 
rious step since ref . does not make any attempt 
to clarify the origin of the extrapolation proce- 
dure and its stability to a change in volume. 

In fact, this size of the lattice was too small to 
yield qualitatively correct results for the quanti- 
ties like < -0V > an d meson masses. These data 
were seriously distorted by the finiteness of the 
volume. The photons, being presumably mass- 
less, can only be more sensitive to such a small 
volume. The same holds for the fermion mass 
since fermions are charged and are, in principle, 
much more sensitive to the presence of the mass- 
less photons then are mesons. Therefore, none 
of the results and the conclusions regarding the 
photon and fermion propagator, can be taken se- 
riously. The quality of the fermion mass and 
renormalized charge data is not good to yield any 
quantitative statements and, if taken seriously, 
they lead to incorrect and inconsistent result. 

For example, consider the fermion mass data 
taken from the DESY group. In general, for 
any mass, Ma, in the scaling region the follow- 
ing equation holds: M A = t"GA(m/t A ). Thus, 
the ratio Mp/M a = G(m/t A ) is a function of 
only one variable. This is true if both masses 
scale. The data for Mp/M a , taken from 0, are 
shown in Fig. 7. Regardless of the value of m, 
the ratio at the critical point, t = 0, always has 
the same value. Thus, the curves R vs (3 have 
to cross at (3 C (insert of Fig. 7). However, data 
show that the curves never cross! The fermion 
mass, as taken from the data, does not scale. It 
doesn't show proper restoration of chiral symme- 
try probably because finite volume effects keeps it 
too heavy. Clearly, if scaling properties were to be 
recovered, the extrapolation would drive us to the 
weaker coupling in order to reduce the fermion 
mass to its value compatible with the phase of re- 
stored symmetry. This is why the conclusions ob- 
tained through extrapolation of the data in ref. [Q 
agreeed so well with renormalized perturbation 
theory - high fermion mass was confused for low 
renormalized charge. This is how the authors of 
ref. misinterpreted their data as supporting the 



zero charge continuum limit. 



Figure 7. Fermion-sigma mass ratio versus bare 
coupling. DESY data. The insert is an illustra- 
tion of the behavior of mass ratio in the scaling 
region. 



As far as the two flavor data is concerned, the 
extensive studies have been done on lattices rang- 
ing from 10 4 to 16 4 . The results have been re- 
ported in ref.|]. The same strategy has been ap- 
plied as in the case of four flavor theory. In Fig.l. 
the mass ratio data are confronted with the pre- 
dictions of the power-law EOS (solid lines). 

A different approach, using microcanonical 
method, is adopted by the Zaragoza group |Q. 
The results for the critical coupling and the crit- 
ical indices obtained by the Illinois and Zaragoza 
group agree. 

A brief comment about perturbation theory 
should be made. The applicability of perturba- 
tion theory in the case of scalar theories has been 
diagnosed early on. For example, in the cr-model 
the weak coupling phase is at low-temperatures. 
Due to the presence of Goldstone bosons, all 
the correlation functions are saturated with the 
massles states and the entire low-temperature 
phase is massless. Every point is a critical point 
in the limit of vanishing magnetic field. Thus, 
the low temperature expansion is an expansion in 
powers of T. The terms of the form exp(— M/T) 
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are absent and there is no danger that they will 
be omitted by using perturbation theory. In this 
way, in principle, the critical region can be ac- 
cessed through perturbation theory which is quite 
accurate close to two dimensions because the crit- 
ical temperature moves to the origin. Clearly 
such reasoning can not be applied to fermions 
simply because the weak coupling phase is sym- 
metric. Thus, no matter how small the cou- 
pling is, perturbation theory omits the Goldstone 
physics as a matter of principle. It can not 
produce bound states that accompany the chiral 
transition and its use is questionable in this case. 

This holds for QED as well although renormal- 
ized charge is bounded from above pJ| . In fact, it 
is difficult to see how perturbation theory could 
communicate with the chiral transition. For ex- 
ample, the fact that the value of the mass ratio 
at the critical point is independent of the bare 
mass is a truely nonperturbative result and could 
not be obtained from any expansion in powers 
of Qffl. Since is dimensionless, at the criti- 
cal point it must be independent of m as well, 
which is clearly impossible in perturbation the- 
ory where screening is sensitive to the fermion 
mass. In other words, conflicting RG trajecto- 
ries must appear if perturbative analysis of the 
renormalized charge is used. 

3. Conclusions 

We have seen so far that the problem of ana- 
lyzing the scaling region of an unknown theory is 
not straightforward. Depending on the accuracy 
of the data, ambiguous results can follow. If (3 C is 
underestimated, and appropriate extrapolations 
done, the "best fit" is given by mean field expo- 
nents supplemented with the log's in the wrong 
place! We have been able to eliminate this pos- 
sibility on theoretical grounds jl2|]. At this stage 
it is fair to say only that the data support the 
power-law scaling with non-gaussian exponents. 
The data of all three groups (|Tof||7|| [{If agree on 
this point. Triviality as realized in 4> theory is 
definitely ruled out. The scenario proposed by 
the DESY group is dismissed on theoretical 
grounds, without a need to consult the data. It 
remains to be seen whether the power law scal- 



ing with 6 ~ 2.3 can be distinguished numeri- 
cally from mean-field plus log's with the log's as 
in Nambu-Jona-Lasinio model. This is the ques- 
tion that needs to be addressed eventually since 
the value of 6 obtained from the simulations and 
is not too far from its mean-field value. Work on 
this problem is in progress. 
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